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Considering the Balitsky-Kovchegov QCD evolution equation in full momentum space, 
we derive the travelling wave solutions expressing the nonlinear saturation constraints 
on the dipole scattering amplitude at non-zero momentum transfer. A phenomenologi- 
cal application to elastic vector meson production shows the compatibility of data with 
the QCD prediction: an enhanced saturation scale at intermediate momentum transfer. 

1 Motivation 

The saturation of parton densities at high energy has been mainly studied for the forward 
dipolc-target scattering amplitude T{r,q = 0, Y), where r, q, Y are, respectively, the dipole 
size, the momentum transfer and the total rapidity of the process. For instance, the corre- 
sponding QCD Balitsky-Kovchegov (BK) equation [2] has been shown to provide a theoret- 
ical insight on the "geometric scaling" properties [3] of the related 7*-proton cross-sections. 
Indeed, it can be related to the existence of a scaling for T(r,q = 0,Y) ~ T(r 2 Q 2 {Y)) 
where the saturation scale is Q 2 (Y) ~ expel" and the constant c can be interpreted as the 
critical speed of "travelling wave" solutions of the nonlinear BK equation [I] . Our theoreti- 
cal and phenomenological subjects are the extension of these properties to the non-forward 
amplitude T(r, q ^ 0, Y), which is phenomenologically relevant for the elastic production of 
vector mesons in deep inelastic scattering. 

2 BK equation in full momentum space 

In order to study the properties of T(r, q 
0, Y), one has first to deal with both concep- 
tual and technical difficulties. It is known 
that the BK formalism has been originally 
derived in impact parameter b but then its 
validity especially at large b is questionable, 
since it leads to non physical power-law 
tails. Hence we start with the formulation 
of the BK equation in momentum q, which 
is more local but has a non-trivial nonlin- 
ear form [5]. In fact, despite this problem, 
the general method of travelling wave solu- 
tions can be extended in the non-forward 
domain [6], It consists in 3 steps: first, one 
solves the equation restricted to its linear 
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part which is related to the non-forward Balitsky Fadin Kuraev Lipatov (BFKL) equation 
[7] for the dipole-dipole amplitude via factorisation and whose solution takes the form of 
a linear superposition of waves. Second, one finds that the nonlinearities act by selecting 
the travelling wave with critical speed c, in a way which, interestingly, is independent of the 
specific structure of the nonlinear damping terms. Third, one obtains after enough rapidity 
evolution, a solution which appears independent from initial conditions (Tq ~ r 270 ) , pro- 
vided these are sharper than the critical travelling wave front profile T ~ r 27c , with 70 > j c - 
Interestingly enough, QCD color transparency satisfies this criterium. Applying these gen- 




Figure 2: p (HI) and <j> (ZEUS) differential cross-sections at W = 75 GeV 

eral results on the non-forward case one finds the following QCD predictions, depending on 
the relative magnitude of three scales involved in the process, namely q, fc^ 1 (the target 
size) and kp 1 = r (the projectile i.e. dipole size). 

• Near-Forward region q <C kp <C kp : Q 2 (Y) ~ k^expcY 

• Intermediate transfer region kp <C q <C kp : Q 2 {Y) ~ q 2 expcY 

• High transfer region q <C kp <C kp : No saturation. 

Our main prediction is thus the validity of the forward travelling wave solution extended 
in the non-forward intermediate-transfer domain but with an enhanced saturation scale by 
the ratio q 2 /k^, where kp is a typically small, nonperturbative scale. Hence we are led to 
predict geometric scaling properties with a purely perturbative initial saturation scale given 
by the transverse momentum. This saturation scale enhancement prediction is confirmed by 
numerical simulations of the BK solutions as shown in Fig[T] 

3 QCD Saturation Model for Exclusive VM production 

The differential cross-section for exclusive vector meson (VM) production at HERA, see 
Fig. 2, can be theoretically obtained from the non-forward dipole-proton amplitude and from 
$J L , the overlap functions between the (longitudinal and transverse) virtual photon and 
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vector meson wave- functions [5] . For completion, we used two different VM wave- functions 
of the literature, without noticeable difference in our conclusions. One writes 

r 1 , 2 

d 2 r / dz & T l(z, r; Q 2 ,M V ) e~ iz ^ T(r, q, Y) , 
Jo 

Following theoretical prescriptions, we consider a forward dipole-proton amplitude Mum 
satisfactorily describing the total DIS cross-sections in a saturation model [9]. We just make 
the saturation scale varying with q 2 , following the trend shown in FigfT] and starting from 
the forward model one Q 2 (Y) 7 one writes 

T{r, q] Y) = 2^R 2 p e- B ^M UM {r 2 Q 2 s (Y,q)) ; Q 2 (q, Y) = Q 2 {Y) (1 + c q 2 ) . 

The factor 2irR 2 e~ Bq comes from the non- 
perturbative proton form factor. For clar- 
ity of the analysis, we considered only B 
and c as free parameters of the non-forward 
paramctrisations, the others being indepen- 
dently fixed by the forward analysis. 

In Table[lJ one displays the x 2 /point ob- 
tained by a fit of p (47 data points) and (j> 
(34 points) total elastic production cross- 
Table 1: Comparison of the x 2 /points sections and of p (50 data points) and <f> (70 

points) differential cross-sections. The Ta- 
ble compares the saturation fits for fixed and g 2 -dependent scales, with a favour for the 
cnhanced-scale model in the total. The model gives a comparable fit with a more conven- 
tional non-saturation model using a Q 2 -dependent slope B oc My + Q 2 - Some of our results 
for the cross-sections are displayed in the figures. In Figj2l one shows the results of the fit 
for p-production (HI) and ^-production (ZEUS) differential cross-sections for a total 7* — p 
energy W = 75GeV and different Q 2 values. Let us finally present our predictions for the 
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Figure 3: Predictions for the DVCS measurements. Left plot: cross-section, right plot: 
differential cross-section. 
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DVCS cross-section, which is obtained without any free parameter from our analysis. In 
FigO they are compared with the available data and the agreement is good in the simple 
chosen parametrisation. 

4 Conclusions 

Let us summarize our new results 

• Saturation at non-zero transfer: The Balitsky-Kovchcgov QCD evolution equation 
involving full momentum transfer predicts (besides the known q = case) saturation in the 
intermediate transfer range, namely for Qq < q < Q, where Qq (resp. Q) is the target (resp. 
projectile) typical scale. 

• Characterisation of the universality class: The universality class of the corresponding 
travelling- wave solutions is governed by a purely perturbative saturation scale Q S {Y) = 
q 2 il(Y), where ft(Y) ~ e cY is the same rapidity evolution factor as in the forward case. 
Consequently the intermediate transfer saturation scale gets enhanced by a factor q 2 /Qq. 

• Phenomenology of Vector mesons: The QCD predictions are applied in the exper- 
imentally accessible intermediate transfer range of vector meson production. The model 
uses an interpolation between the forward and non-forward saturation scale together with a 
parameter-frozen forward saturation model. It fits better the data on p (HI) and cj> (ZEUS) 
cross-sections than for a non-enhanced saturation. 

• Prospects: The next phenomenological prospect is to add charm to the discussion, 
both with the modification of the forward case by including the charm contribution [lOj 
and by also considering the production of \P mesons. On a theoretical ground, it would be 
interesting to go beyond the mean-field approximation of the BK equation. 
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